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"Remote Sensing Reflectance and Inherent Optical Properties in the Chesapealce Bay " 


Abstract: 

We used an extensive set of bio-optical data and radiative transfer (RT) model 
simulations of radiation fields to investigate relationships between inherent optical properties and 
remotely sensed quantities in the optically complex, mid-mesohaline Chesapeake Bay waters. 
Field observations showed that the chlorophyll algorithms used by the MODIS (MODerate 
resolution Imaging Spectroradiometer) ocean color sensor (i.e. Chlor_a, chlor_MODIS, 
chlor_a_3 products) do not perform accurately in these Case 2 waters. This is because, when 
applied to waters with high concentrations of chlorophyll, all MODIS algorithms are based on 
empirical relationships between chlorophyll concentration and blue-green wavelength remote 
sensing reflectance (Rrs) ratios that do not account for the typically strong blue-wavelength 
absorption by non-covarying, dissolved and non-algal particulate components. Stronger 
correlation was observed between chlorophyll concentration and Rrs ratios in the red (i.e. 
i?rs(677)/j?rs(554)) where dissolved and non-algal particulate absorption become exponentially 
smaller. Regionally-specific algorithms that are based on the phytoplankton optical properties in 
the red wavelength region provide a better basis for satellite monitoring of phytoplankton blooms 
in these Case 2 waters. Good optical closure was obtained between independently measured Rrs 
spectra and the optical properties of backscattering, and absorption, a, over the wide range of 
in- water conditions observed in the Chesapeake Bay. Observed variability in the quantity f/Q 
(proportionality factor in the relationship between Rrs and the water inherent optical properties 
ratio b\>l(a+b\,)) was consistent with RT model calculations for the specific measurement 
geometry and water bio-optical characteristics. Data and model results showed that f/Q values in 
these Case 2 coastal waters are not considerably different from those estimated in previous 
studies for Case 1 waters. Variation in surface backscattering significantly affected Rrs magnitude 
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across the visible spectrum and was most strongly correlated (R 2 =0.88) with observed variability 
in Ris at 670 nm. Surface values of particulate backscattering were strongly correlated with non- 
algal particulate absorption, a nap , in the blue wavelengths (R 2 =0.83). These results, along with 
the measured values of backscattering fraction magnitude and non-algal particulate absorption 
spectral slope, suggest that suspended non-algal particles with high inorganic content are the 
major water constituents regulating b\, variability in the mid-mesohaline Chesapeake Bay. 
Remote retrieval of surface bh and a nap from #rs(670) can be used in regionally-specific satellite 
algorithms to separate contribution by non-algal particles and dissolved organic matter to total 
light absorption in the blue, and monitor non-algal suspended particle concentration and 
distribution in these Case 2 waters. 

Keywords: bio-optics; satellite retrievals; radiative transfer; Case 2; estuaries; coastal waters 
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1. Introduction 

Changes in the concentrations and distribution of organic and inorganic, particulate and 
dissolved substances are of major water quality and ecologic concern in the Chesapeake Bay 
estuary. Increased nutrient loadings in recent years, driven mainly by human population growth 
and changes in land use, have increased phytoplankton concentrations above historical levels 
(Harding and Perry, 1997). Changes in particulate and dissolved substances are linked to 
important processes in the estuary, such as freshwater river discharges, nutrient and light 
availability in the water column, tidal mixing, bottom resuspension, and microbial activity (e.g. 
Kemp and Boynton, 1992; Malone, 1992). Because composition and concentrations of seawater 
constituents influence water optical characteristics, in-situ measurements of optical properties 
have been made by several ship-based monitoring programs during the last decades to examine 
water quality and assess progress in reversing eutrophication in the Bay (e.g. Glibert et al., 1995; 
Johnson et al., 2001). By using appropriate bio-optical algorithms, remote sensing of “ocean 
color” offers the capability of extending field observations beyond the restricted in-situ sampling 
domain. Several aircraft ocean-color instruments (using both “passive” and “active” systems) 
have been used to remotely measure surface chlorophyll-a concentrations [chl-a] and determine 
changes in phytoplankton biomass in the Chesapeake Bay (e.g. Hoge and Swift, 1981; Harding 
et al., 1992; Lobitz et al., 1998). However, only a limited number of studies has been published 
on the interpretation of satellite ocean color imagery and the applicability of currently used 
satellite algorithms for these optically complex estuarine waters (e.g. Harding et al., 2005). 

In order to use satellite ocean color to extract information on water composition it is 
necessary to develop effective bio-optical algorithms relating the remotely sensed water 
reflectance either directly to surface concentrations of optically active water constituents 
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(empirical algorithms; e.g. Clark, 1997; O’Reilly et al., 2000), or to their optical properties based 
on principles derived from radiative transfer theory (semi-analytical inversion models; e.g. 
Garver and Siegel, 1997; Maritorena et al., 2002). Since the launch of the Coastal Zone Color 
Scanner (CZCS) in October 1978, satellite ocean color has contributed significantly to gaining a 
better understanding of biological activity in open-ocean “Case 1” waters where phytoplankton 
and co-varying material are the major optical components (e.g. Yentsch, 1993; Longhurst et ah, 
1995; Gregg and Conkright, 2001). With more channels in the visible part of the spectrum, 
CZCS’s follow-on sensor SeaWiFS (Sea-viewing Wide Field of view Sensor; launched in 1997) 
and the newer instrument MODIS (MODerate resolution Imaging Spectroradiometer, launched 
in 1999) allowed for further improvements in satellite retrievals of biogeochemical variables in 
Case 1 waters (e.g. Yoder and Kennelly, 2003; Carder et ah, 2004; McClain et ah, 2004). A more 
difficult challenge, however, has been developing bio-optical algorithms suitable for use in 
optically complex “Case 2” waters, like the Chesapeake Bay, where multiple, co-existing but not 
necessarily co-varying, dissolved and particulate, marine and terrigenous substances affect ocean 
color (e.g. Morel and Prieur, 1977; Carder et ah, 1991 ; Darecki and Stramski, 2004). 

Harding et ah (2005) used SeaWiFS [chl-a] data and in-situ measurements of water 
reflectance and [chl-a] to examine the accuracy of satellite observations and the applicability of 
the SeaWiFS empirical chlorophyll algorithm OC4v4 (O’Reilly et ah, 2000) for the Chesapeake 
Bay and Middle Atlantic Bight. They found that SeaWiFS reliably captured seasonal and inter- 
annual variability of phytoplankton biomass in the lower Bay. However, the OC4v4 algorithm 
significantly overestimated chlorophyll in the mid-mesohaline and upper Chesapeake Bay, due to 
strong absorption by dissolved organic matter and non-algal particles that are not sufficiently 
accounted for in this empirical chlorophyll algorithm (Harding et ah, 2005). 
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These studies underscore the need to develop more accurate Case 2 algorithms based on 
detailed in-situ bio-optical characteristics. However, certain optical properties of the Chesapeake 
Bay remain poorly characterized. Magnuson et al. (2004) used an extensive set of bio-optical 
data to examine the parameterization and validation of the semi-analytical Garver-Siegel- 
Maritorena (GSM01) model (Maritorena et al., 2002) for application in the Chesapeake Bay and 
the Mid- Atlantic Bight. A major limitation in their efforts was that information on the 

backscattering (Z>b) characteristics of Chesapeake Bay is extremely scarce. Due to the lack of in- 

•x 

situ data, modeling of b\, is often based on assumptions regarding scattering angular shape, 
backscattering magnitude and spectral dependence. According to Magnuson et al. (2004) the lack 
of b\, measurements in the Bay affected backscattering parameterizations in their study and posed 
a significant limitation on the evaluation of the GSM01 model’s backscattering product. 
Tzortziou et al. (submitted) recently applied detailed in-situ data and radiative transfer 
simulations to examine bb variability and modeling of backscattering processes for the mid- 
mesohaline Chesapeake Bay. Information is scarce on the contribution of various water 
constituents to the backscattering characteristics of Chesapeake Bay, and the application of 
remotely sensed ocean color data to the remote retrieval of particulate backscattering variability 
in these Case 2 waters. 

Previous studies have shown that the remote sensing reflectance Rrs, defined as the 
upwelling radiance emerging from the ocean divided by the downwelling irradiance reaching the 
water surface, is, to a first approximation, proportional to the ratio of water backscattering, b b , 
and absorption, a, coefficients, bb/(a+bb) (e.g. Morel and Prieur, 1977). However, the magnitude 
and geometrical structure (or shape) of the upwelling radiance field within the ocean depend on 
both the water inherent optical properties (IOPs; e.g. volume scattering function (VSF), bb, a), as 
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well as on the processes and environmental variables (e.g. solar zenith angle, aerosol load, 
surface waves) that affect illumination conditions (Austin, 1974; Morel and Gentili, 1991; 1993). 
Measurements, theoretical studies, and radiative transfer computations of the effects of these 
parameters on the relationship between Rrs and the IOPs ratio b b /(a+b b ) have been performed for 
Case 1 waters (e.g. Voss, 1989; Morel and Mueller, 2002; Voss and Morel, 2005). Further 
studies are needed for Case 2 waters, like the Chesapeake Bay, that are characterized by high 
bio-optical complexity. 

In this paper, we use radiative transfer (RT) modeling and detailed in-situ bio-optical 
measurements, including new data on particulate backscattering, to examine relationships 
between remotely sensed Rrs spectra and IOPs for improved remote retrievals of biogeochemical 
variables in the mid-mesohaline Chesapeake Bay. Our main objectives were to: (i) assess the 
applicability of the chlorophyll algorithms currently used by SeaWiFS’ successor MODIS, and 
examine alternative chlorophyll algorithms based on the red wavelengths where absorption by 
non-algal particles and dissolved matter are minimal; (ii) evaluate variability in the relationship 
between Rrs and b\J(a+b b ) for our measurement geometry, and compare our observations with 
RT calculations for Case 2 waters and previous studies for Case 1 waters; (iii) determine the 
relative roles of different particulate substances in light backscattering in the mid-mesohaline 
Chesapeake Bay, and examine methods for remote retrieval of surface b b and suspended particles 
in these waters. 

2. Theoretical background 

Remote sensing of ocean color relies on detecting the light signal that leaves the water 
surface and reaches a sensor onboard a satellite, carrying with it information on water IOPs. The 
magnitude of the ocean remote sensing reflectance Rrs is related to the processes of 
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backscattering, which redirects downwelling photons to travel upward and eventually leave the 
water surface, and absorption, which converts photons energy to other forms of energy such as 
heat or chemical energy. As upward traveling photons exit the water they interact with the air- 
water interface by refraction and internal reflection. Rrs can, therefore, be related to 
backscattering and absorption according to (e.g. Austin, 1974; Preisendorfer, 1976; Gordon et 
ah, 1975; 1988; Lee etal., 1994): 


Rrs(X) = ■ - ( - ’ a - ) V ) - ^ 

Q n w 2 a(A) + b b (A) 


0 ) 


where, bb(X) is the total backscattering coefficient at wavelength X, a(X) is the total absorption 
coefficient, /( W>a ) is the water-air transmittance, t( a, W ) is the transmittance from air to water, and n w 
is the real part of the water refractive index. The quantity /is a complex function of wavelength, 
water IOPs (single scattering albedo and volume scattering function), solar zenith angle ( 0 O ), 
aerosol optical thickness, and wind speed (Gordon et al., 1975; Kirk, 1984; Morel and Mueller, 
2002). The quantity Q is defined as the ratio of upwelling irradiance to upwelling radiance, 

Q= Eu(X)/Lu(X) (Austin, 1974). Therefore, Q is influenced by the processes and environmental 
variables that affect the geometrical structure of the anisotropic upward radiance field (i.e. IOPs 
and environmental variables mentioned above, direction of upward traveling photons, and 
azimuth angle) (Morel and Mueller, 2002). 

Variability in the ratio jlQ and its implications in ocean color remote sensing have been 
the focus of several theoretical and RT computational studies. Early studies showed that the 
quantity /has an average value of about 0.32-0.33 (Gordon et al., 1975; Morel and Prieur, 1977) 
when the sun is near zenith. However, according to Morel and Mueller (2002) the global range of 
variation in /is from about 0.3 to 0.6. Morel and Gentili (1993) found that Q values generally 
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range from 3 to 6. According to their study, both the / and Q functions experience concomitant 
increases when the sun zenith angle increases. Therefore, the jlQ ratio would be expected to be 
less dependent on solar zenith angle than either f ox Q alone. Studies by Morel and Gentili (1993) 
and Gordon et al., (1988) have shown that the ratio f/Q is relatively independent of solar zenith 
angle for sun and viewing angles expected for the MODIS orbit, with average f/Q = 0.0936, 
0.0944, 0.0929, and 0.0881, (SD = 0.005), for X= 440, 500, 565, and 665 nm, respectively 
according to Morel and Gentili (1993), and f/Q = 0.0949, for 6 0 < 20° according to Gordon et al. 
(1988). Thus, in various bio-optical algorithms that relate Rxs to \Q?s,f/Q is assumed to be 
independent of wavelength and solar zenith angle (e.g. the MODIS semi-analytical chlorophyll 
algorithm; Carder et al., 2002). 

Using detailed RT computations Morel and Mueller (2002) showed that f/Q varies within 
the range 0.08 - 0.15. For simplicity. Morel and Mueller performed their RT calculations for 
homogeneous Case 1 waters, which allowed them to model all water IOPs (a, 6b, and VSF) as a 
function of the chlorophyll-a concentration (Morel et al., 2002). Their theoretical computations 
of Q showed very good agreement with measurements of radiance distribution performed by 
Voss and Morel (2005) over a large range of chlorophyll concentrations in the Case 1 waters of 
Baja California. The degree of agreement between theoretical computations and measurements 
of f/Q variability in Case 2 waters, however, remains unknown. One of the main objectives in 
this paper is to examine variability in the quantity f/Q and, as a result, in the relationship of 
Rxs : 6b/(a+6b), for the geometry of our measurements, and over a wide range of bio-optical and 
environmental conditions for the Case 2 Chesapeake Bay waters. 
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3. Methods 

3.1. Location and duration of measurements 

Measurements of water inherent optical properties and radiation fields were performed at 
four stations designated: Poplar Island (PI), Herring Bay (HB), Tilghman Island (TI), and Jetta 
(JT), in the mid-mesohaline Chesapeake Bay estuary, on the eastern coast of the US (38.71° - 
38.89° N latitude, 76.34° - 76.54° W longitude) (Figure 1). Measurements were made during 17 
cruises performed in the spring, summer, and fall months between June 2001 and November 
2002 (Table 1). 

3. 2. In-situ and laboratory measurements of water optical characteristics 

In-situ vertical profiles of total (minus pure water) attenuation, c t . w (X, z), and absorption, 
a,. w (k, z), were measured at nine visible wavelengths using a WETLabs Spectral AC-9 (Table 2). 
Temperature and salinity were measured with Hydrolab Datasonde 4a, and these data were used 
to correct AC9 measurements for the temperature and salinity dependence of absorption by pure 
water (Pegau et al., 1997). Total particulate scattering, b p (k, z), was estimated as the difference 
between c,. w (k, z) and a t . w (k, z), after applying additional corrections to account for scattering 
losses manifested as overestimates of measured absorption (Kirk, 1992; Tzortziou et al., 2004). 
An ECOVSF3 instrument (WetLabs Inc; Moore et al., 2000) was used to measure backscattering 
at three angles (100 °, 125 °, 150°) and three visible wavelengths (450, 530, 650 nm). 
Measurements were corrected for attenuation and were integrated (90-180°) to obtain the total 
backscattering coefficient according to manufacturer’s instructions. Measurements by Boss et al. 
(2004) in the Case 2 waters off the coast of New Jersey showed that estimates of b b using an 
ECOVSF3 instrument were in very good agreement (R 2 =0.99) with b\, measurements using a 
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HOBILabs Hydroscat-6 (Maffione and Dana, 1997). These results increase confidence on the 
accuracy of the backscattering measurements technique, especially since the instruments have 
large differences in design and calibration (Boss et al. 2004). 

Water samples collected from discrete depths were filtered using 0.22 pm pore-diameter 
polycarbonate membrane filters (Poretics) to separate the dissolved material. Suspended 
particulate matter was collected on 25 mm glass fiber filters (Whatman GF/F). Absorbance 
spectra of chromophoric dissolved organic matter (CDOM), phytoplankton, and non-algal 
particulate matter were measured using a Cary-IV dual-beam spectrophotometer using the 
approach of Gallegos and Neale (2002). Particulate absorbance spectra were corrected for 
scattering errors using a pathlength amplification factor estimated empirically by comparing 
particulate optical density measured on filters and in particle suspension (Mitchell et al., 2000; 
Tzortziou, 2004). Absorption spectral slope coefficients, describing the exponential decline of 
absorption with increasing wavelength for CDOM, Scdom> and non-algal particulate matter, S nap , 
were determined by applying non-linear least-squares regression to the absorption coefficients 
versus wavelength (290-750 nm) (Blough and Del Vecchio, 2002). Chlorophyll-a concentrations 
were measured spectrophotometrically in 90% acetone extracts of seston collected on 47 mm 
GF/F glass fiber filters (Jeffrey and Humphrey, 1975). 

3.3. In-situ data and theoretical estimations of radiation fields 

Two sensor arrays were used to measure underwater radiation fields in this work, 
depending on instrument availability. On eight cmises (Table 1) underwater upwelling, £u(z), 
and downwelling, £d(z), spectral irradiance profiles, and above-water surface downwelling 
irradiance, Es, were measured using Satlantic 00-200 seven-channel irradiance sensors (Table 
2). The irradiance sensors were mounted on a custom frame so that up- and down-welling 
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sensors were almost coplanar. On the rest of our cruises we used a Satlantic MicroPro free- 
falling radiometer (Table 1) to measure water column profiles of upwelling radiance, Zu(z), and 
Ed(z) in the water column, and Satlantic OCR-507 Surface Reference Irradiance sensors for 
measurements of Es, at 14 wavebands (Table 2) . The optical heads of MicroPro have relatively 
small diameters (6.4 cm) compared to other radiometric instruments and are less subject to 
instrument self-shading (Harding and Magnuson, 2002). The instrument contains a pressure 
sensor that gives depth, and a miniature biaxial clinometer for tilt measurements (accuracy of 
0.2°) (Satlantic MicroPro operation manual, 2002). 

When using the Micropro, three casts were made at each station, and all casts were 
completed within 5-8 minutes. Measurements that did not meet quality control (e.g. casts 
characterized by large tilt-angles or changing cloudiness conditions) were omitted from analysis. 
In cases when more than one cast were of good quality Z,u was estimated as the average of all 
casts. Measurements of Zu(z) were corrected for instrument self-shading errors (Gordon and 
Ding, 1992; Zibordi and Ferrari, 1995) and for the depth offset between the Ed and Lu sensors. A 
correction was applied to the radiometric measurements through the instrument’s calibration for 
the immersion effect (Satlantic MicroPro operation manual, 2002). 

To estimate the water leaving radiance Zw, we extrapolated underwater Zu(z) 
measurements to just below the water surface z=0‘ and estimated transmittance through the 
water-air interface. The upwelling radiance just below the water surface, Zu(0',X) was estimated 
through non-linear least squares fitting (SigmaStat software) of measured Zu(z,Z) to the equation 
Zu(z,Z)= Zu(0',Z) exp(-Ki u • z ), where Ki u is the diffuse attenuation coefficient for the upwelling 
radiance, to a good approximation constant in the layer 0-3 m used in the regressions. The 
coefficients of determination (R 2 values) for the non-linear exponential fits were in most cases 
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better than 0.99. To estimate Zw, we calculated the propagation of Zu(0~,Z) through the water-air 
interface: 

Lw{X,0,(p) = Lu(O-,JL,0',p) Q- r ~%’?9 (2) 

where 6’ is the direction of the upward traveling photons incident from the water body onto the 
water surface, 9\s the direction of the transmitted photons, r(9’,6) is the Fresnel reflectance for 
the associated directions 9‘ and 9, and n w is the index of refraction of water (« w « 1.34) (Mobley, 

1 994). For the geometry of our measurements, the zenith angle of water-leaving radiance and the 
nadir angle of in- water upward radiance are zero (9-9=0) and the transmittance is 
(l-r(i9’,#))«0.98 (Mobley, 1994). Therefore, Zw(Z) can be estimated from Eq. 2, as: 

Zw (Z) = 0.544 Zu(0\ X) (3) 

Based on its definition, Rrs(X) was estimated as the ratio of the water-leaving radiance Zw(Z) and 
measured surface downwelling irradiance Es{X). Eq. 1 and measured Rrs and IOPs (a and 6b) 
were used to estimate variability in the quantity f/Q. 

For those cases when in-situ measurements of Zu profiles were not available we used the 
extensively validated Hydrolight underwater radiative transfer (RT) program (Mobley, 1988; 
1994) to estimate Zu(O'), Zw, and Rrs for the bio-optical conditions observed in the Bay. 
Measured IOPs (i.e. a, 6, c, 6b) and above-water Es were used as inputs to perform the model 
calculations. The Pope and Fry (1997) absorption values for pure water and the Smith and Baker 
(1981) seawater scattering coefficients with the Rayleigh-like pure- water scattering phase 
function were used in our RT simulations. Raman scattering, and CDOM and chlorophyll-a 
fluorescence were included in all model runs. Tzortziou et al. (submitted) discuss details of the 
methodology and results of RT model simulations using Hydrolight for the mid-mesohaline 
Chesapeake Bay, and the good RT optical closure they obtained. 
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3.4. MODIS chlorophyll algorithms and products 

Three algorithms have been used with MODIS data for remote retrieval of chlorophyll 
concentrations: 

i) The SeaWiFS-analog OC3M chlorophyll algorithm (currently operational MODIS level-2 
standard chlorophyll product ‘Chlor-a’) is an empirical algorithm that relates the maximum of 
three blue-green Rxs ratios to [chl-a] (O’Reilly et al., 2000; Eq. 4 in Table 3). 

ii) The MODIS empirical-HPLC chlorophyll algorithm (‘chlor_MODIS’ product), is based on an 
empirical relationship derived from high performance liquid chromatography (HPLC) 
measurements of [chl-a] and blue-green Rxs ratios measured in Case 1 and Case 2 waters (Clark, 
1997; updated 19 Feb. 2002, D. Clark, personal communication) (Eq. 5 in Table 3). A recently 
proposed version of the algorithm (updated 19 March 2003, D. Clark, personal communication) 
is a 5 th - order polynomial, expected to perform better in very high and very low chlorophyll 
environments. 

iii) The MODIS semi-analytical chlorophyll algorithm (‘chlor_a_3’ product) is based on a bio- 
optical model that relates Rxs to backscattering and absorption by phytoplankton and gelbstoff 
(combined term for CDOM and non-algal particles). Due to the similarity in the absorption 
spectral shapes of CDOM and non-algal particles the algorithm cannot separate contribution of 
each component to ‘gelbstoff absorption (Carder et al., 2002). The algorithm uses a semi- 
analytical approach in the sense that the relationship between i?rs and water IOPs is based on 
radiative transfer theory and the approximate formula Eq. 1 . However, for tf P hyt(675) larger than 
0.03 m' 1 , or [chl-a] higher than 1. 5-2.0 mg m' 3 , the MODIS semi-analytical algorithm 
automatically switches to an empirical relationship between [chl-a] and blue-green Rxs ratios 
(Eq. 6 in Table 3). 
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In summary, when applied to waters with [chl-a] larger than ~2 mg m" 3 > which is 
typically the case in coastal waters like Chesapeake Bay, all three MODIS algorithms are based 
on empirical relationships between [chl-a] and Rrs ratios at blue-green wavelengths. To examine 
the performance of the MODIS algorithms in the mid-mesohaline Chesapeake Bay region we 
applied the satellite algorithms (Eq. 4, 5, and 6 in Table 3) to Rrs spectra estimated based on our 
in-situ bio-optical data. We then compared these [chl-a] estimates to measured surface 
chlorophyll concentrations. 
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4. Results 

4.1. Water optical characteristics and RT model simulations 

Water bio-optical characteristics showed large spatial and temporal variability during our 
measurements in the mid-mesohaline Chesapeake Bay. Relatively clear waters were observed 
during the fall cruises, with low biological activity and relatively low total absorption, total 
attenuation and chlorophyll concentration values. Higher chlorophyll concentrations, associated 
with large surface phytoplankton bloom events, were observed in spring and summer (Tzortziou, 
2004). Surface [chl-a] values ranged from 3.5 mg m' 3 (station PI, 13 November 2001) to 
74 mg m' 3 (station HB, 1 1 June 2001) with an average value of 14.7 mg m' 3 (Figure 2). 

Measurements of absorption properties showed that non-covarying CDOM and non-algal 
particles in these Case 2 waters contribute considerably to light attenuation at the blue 
wavelengths used in the MODIS chlorophyll retrievals. Combined contribution by CDOM and 
non-algal particles to surface (0-lm) a t -w was on average 59% at 488 nm. Absorption by non- 
algal particles alone was on average 41% of a t . w at 488 nm, getting as large as 56% of a t . w (488) 
at the highly turbid station JT on 9 July 2001. Contribution by CDOM and non-algal particles to 
surface a t . w was even larger at shorter wavelengths (e.g. 73% at 412 nm) as a result of the 
exponential increase in the absorption of both substances with decreasing wavelength (Bricaud et 
al., 1981; Kishino et al., 1985). Estimated Scdom had an average value of 0.018 nm' 1 with a 
standard deviation of 0.0032 nm" 1 . S nap varied in a relatively narrower range compared to Scdom, 
and had an average value of 0.01 1 nm' 1 and a standard deviation of only 0.001 nm' 1 . 

Large spatial and temporal variation was observed in the magnitude of particulate 
backscattering in the Bay, with surface &b(530) ranging between 0.013 - 0.166 m' 1 (Tzortziou, 
2004). Higher b\> values were observed consistently at the near shore site JT, compared to the 
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other three stations. Estimated backscattering fraction values, bfb, ranged from 0.006 to 0.036 at 
530 nm, with the largest b\Jb value of 0.036 occurring close to the bottom at station JT. Average 
bjb values, estimated as the average of measurements at all depths at the four stations, were 
0.0133 (SD = 0.0032), 0.0128 (SD = 0.0032), and 0.0106 (SD = 0.0029), at 450, 530 and 650 
nm, respectively. Surface bjb at 530 nm, estimated from measurements within 1 m from the 
surface, had an average value of 0.0125 (SD = 0.003) at the four stations. Surface &t>/&(530) was 
typically higher at the turbid station JT, with an average value of 0.0146 (SD = 0.0037). 

Observed variability in water IOPs resulted in large spatial and temporal variability in 
the magnitude of measured 7?rs spectra (Tzortziou et al., 2004). However, in all cases, maximum 
values of i?rs occurred at the green wavelengths because of the large pure-water absorption in the 
red, and the large CDOM and non-algal particulate absorption in the blue region of the spectrum. 

4.2. MOD1S chlorophyll-a algorithm performance 

We examined the performance of the MODIS chlorophyll algorithms for the mid- 
mesohaline Chesapeake Bay by applying the MODIS algorithms to Rr s spectra derived from our 
in-situ bio-optical data, and comparing estimated [chl-a] products to our in-situ measurements of 
surface chlorophyll concentrations. For the water conditions we observed in the mid-mesohaline 
Chesapeake Bay chlorophyll concentrations were consistently higher than 3.5 mg m' . For these 
[chl-a] values all MODIS chlorophyll retrievals are based on empirical relationships between 
[chl-a] and blue-green Rts ratios (Eq. 4-6, Table 3). Comparisons between estimated [chl-a] 
values with in-situ measured [chl-a] showed large disagreements for all satellite algorithms 
(Figure 3), although some of the algorithms were expected to perform better in high chlorophyll 
environments (e.g. 5 th order polynomial, “chlor_MODIS” product). Coefficients of determination 
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were in all cases between 0.3 and 0.35 (Table 4). These results were in agreement with the 
relatively weak correlation we observed between measured surface logio([chl-a]) and the Rxs 
ratios used in the MODIS chlorophyll retrievals (i.e. logio(i?rs(443)//?rs(554)) and 
logio(/?rs(488)/i?rs(554)). Coefficients of determination in the applied linear (and cubic) 
regressions were less than 0.39 (Table 4). 

We then examined correlation between surface [chl-a] and various two-band Rxs ratios 
using MODIS wavebands in the blue, green, and red wavelength regions (Table 4). The strongest 
correlation was found between [chl-a] and the ratio i?rs(677)/i?rs(554) (R 2 =0.54). The empirical 
relationship found between [chl-a] and the red-green band ratio was consistent with expectations 
based on our in-situ bio-optical data (Tzortziou, 2004). Based on Eq. 1, and assuming 
fl Q' (t-w-a ta-w) /n w 2 independent of wavelength and b b « a (Carder et al., 2002), the Rxs ratio 
677:554 is proportional to the backscattering ratio and inversely proportional to the absorption 
ratio at these two wavelengths (i.e. i?rs(677)/i?rs(554) = Z>b(677)/Z>b(554) • a(554)/a(677)). For the 
mid-mesohaline Chesapeake Bay, b\, in the green was strongly correlated with 6 b in the red 
(R =0.99; Table 4). Strong correlation and linear relationship (R =0.92) was also found between 
a t . w (554) and a t - w (677) based on our AC9 data (Table 4). As phytoplankton is the major absorber 
in the 677 nm wavelength region (other than pure water), absorption at 677 was expected to be 
strongly correlated with [chl-a], with some variation due to changes in phytoplankton species 
composition, physiological state, and size (Bricaud et al., 1995). Indeed, least squares regression 
for our data gave: a t . w (677) = 0.0166 • [chl-a] + 0.0603 (R 2 = 0.92), with the small intercept due 
to the small, residual absorption by CDOM and non-algal particulate matter at 677 nm. Including 
the effect of chlorophyll fluorescence at 677 nm using detailed RT model simulations (Tzortziou, 
2004), we found that the relationship observed between [chl-a] and 7!rs(677)/i?rs(554) was 
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consistent with predictions based on independently derived relationships between [chl-a] and the 
ratio of water IOPs (i.e. b b , a) at the two specific wavelengths. 

4.3. Observed variability in f/Q 

We examined the relationship between Rrs{\) and the ratio of water IOPs ^0) j n 

a(/l )+b b {X) 

the mid- mesohaline Chesapeake Bay using measured and model-estimated Rrs spectra and in- 
situ data of surface b b (X) and a(X) to: (i) examine the closure between independently measured 
inherent and apparent optical properties based on Eq. 1 ; (ii) assess variability in the quantity f/Q 
(from the slope between Rrs and b b /(a+b b )) for the nadir- vie wing geometry of our measurements; 
and (iii) compare field results with RT model simulations for these Case 2 waters and previous 
theoretical studies on f/Q variability for same viewing geometry but Case 1 waters (Morel and 
Mueller, 2002). 

Rrs values at 443, 532 and 670 nm were strongly correlated with b b /(a+b b ) at all three 
wavelengths, whether Rrs was measured by the Satlantic MicroPro (Figure 4, solid circles) or 
calculated by Hydrolight when MicroPro was not available (Figure 4, open circles). Coefficients 
of determination (R 2 values) were larger than 0.9 in all cases. The linear least-square regression 
fits for Rrs measured versus b\J(a+b\y ) were very similar to the linear least-square regression fits for 
Rrs model versus bf(a+b\/), suggesting a good agreement between data and RT model simulations. 
ft t 

1 — ( "' ,a) (a ’“ ,) , estimated as the slope coefficient in the linear least square regression of f?rs me asured 

Q nf 

or i?rs m0 dei versus bf(a+b b ) (Eq. 1), ranged from 0.0485 to 0.058 for different wavelengths with 
coefficients of variation in the estimated slopes less than 2.5%. Estimates of the quantity 
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f t t 

— — — { *- w) derived as the ratio Rrs : b b /(a+b b ) for each individual data set of Rrs and IOPs had 
Q n K 

a standard deviation of less than 10% at all wavelengths, showing that the modeled f/Q (i.e. slope 
in the regression Rrs vs b b /(a+b b )) represented well the individual measurements over the whole 
range of conditions. 

4.4. Relationships between bt,, particulate matter, and Rrs 

To examine the effect of particulate backscattering on Rrs variability for the mid- 
mesohaline Chesapeake Bay we applied simple linear regression to the Rrs and b b values 
measured over the observed wide range of in-water optical characteristics (Figures 5(a)-(c)). 

Both measured and model estimated Rrs values at 670 nm were strongly correlated with surface 
hb at 650 nm (Figure 5(c)) (R 2 =0.88). The correlation between Rrs and b b was statistically 
significant but not as strong at the shorter wavelengths (i.e. i?rs(443) versus 6b(450) in Figure 
5(a) with R 2 =0.67, and i?rs(532) versus hb(530) in Figure 5(b) with R 2 =0.77). 

Variability in i?rs(670) in these waters is mainly driven by changes in particulate 
backscattering because of the relatively large contribution by pure water to total absorption at 
670 nm. Surface b b in the red ranged between 0.008 and 0.12 m' 1 during our measurements in the 
Bay. This variability in backscattering, by more than an order of magnitude, corresponds to more 
than an order of magnitude changes in i?rs(670) (based on Eq. 1). Rrs is also affected by changes 
in total in-water absorption. However, changes in dissolved and particulate absorption at 670 nm 
by an order of magnitude (i.e. a t . w (670) in the range 0.1-1 m' 1 for the mid-mesohaline Bay) are 
masked by the relatively large pure water absorption (a w (670)= 0.44 m' 1 ; Pope and Fry, 1997) 
translating to changes in total absorption by less than a factor of 3. At shorter visible 
wavelengths, interference by pure water absorption is minimal. Therefore, Rrs at 443 and 530 nm 
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are affected strongly by changes in both total absorption (strong contribution by non-covarying 
particulate and dissolved components) and total backscattering (contribution by suspended 
particles). 

Surface particulate backscattering at all three wavelengths 440, 530 and 650 nm (only 
Z> b (650) nm shown in Figure 6) was strongly correlated with absorption by non-algal particles in 
the blue (e.g. 412 nm) and UV (e.g. 380 nm) where non-algal particulate absorption is 
particularly strong (R 2 = 0.83 for the linear least-square regression bb(650) vs a nap (380); Figure 
6). The relationship between surface [chl-a] and surface 6 b was highly variable (Figure 7) and 
correlation was considerably smaller (R 2 =0.57) compared to that between 6 b and a nap . 

Given the above results, a strong correlation would be expected between i?rs(670) and 
magnitude of absorption by non-algal particles. Indeed, coefficients of determination for the 
linear least-squares regression between i?rs(670) and surface a nap were 0.7 and 0.74 for 
absorption measurements at 412 and 380 nm, respectively (Figures 8(a)-(b), Table 4). 
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5. Discussion 

Obtaining accurate information on the composition and concentration of dissolved and 
particulate, organic and inorganic matter in near shore waters using remotely sensed ocean color 
is critical for primary production studies, coastal water quality monitoring, and carbon cycling 
modeling. However, application of currently operational satellite algorithms in near shore waters 
like the Chesapeake Bay often results in erroneous retrievals, as shown here for the MODIS 
chlorophyll algorithms. In this paper we applied an extensive bio-optical dataset and detailed RT 
model results to examine the relative roles of different substances in light absorption and 
backscattering in the mid-mesohaline Chesapeake Bay, derive relationships between Rrs and 
optical properties of individual water components, and evaluate methods for remote retrieval of 
surface [chl-a], b\» and suspended particles in these optically complex Case 2 waters. 

Effective remote retrieval of biogeochemical variables in coastal waters depends to a 
large extent on the accuracy of, and consistency among, the in-situ data used in the development, 
validation, and improvement of the applied remote sensing bio-optical algorithms. Thus, testing 
the degree of closure among bio-optical quantities independently measured in the field becomes 
critical for remote sensing applications. Before applying our bio-optical data to the interpretation 
of satellite ocean-color observations and the investigation of new bio-optical relationships for the 
Chesapeake Bay, we examined the consistency and degree of optical closure among our 
measurements using RT modeling (Tzortziou et al., 2004). Close agreement between data and 
model results over a wide range of conditions, with average absolute differences between 
modeled and measured Rxs less than 1 0%, demonstrated very good optical closure between 
independently measured quantities. These results increased confidence in the accuracy of, and 
consistency among, our in-situ bio-optical data. 
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5.1. MODIS chlorophyll-a algorithm performance 

Chlorophyll concentrations larger than 2 mg m" 3 are typically observed in the mid- 
mesohaline Chesapeake Bay (e.g. Magnuson et al., 2004). During our measurements, which 
covered a wide range of in-water conditions, measured [chl-a] was consistently larger than 3.5 
mg m' . For such conditions ([chl-a]>2 mg m ) all MODIS chlorophyll retrievals are based on 
empirical relationships between [chl-a] and blue-green reflectance ratios (i.e. 443, 488 nm 
MODIS bands). When we applied these algorithms to measured Rrs spectra to examine the 
MODIS algorithm performance in the mid-mesohaline Bay, we found large disagreement 
between estimated and measured surface [chl-a] (Figure 3). The weak performance of MODIS 
chlorophyll retrievals in these Case 2 waters is not surprising considering our results showing 
large contribution by non-covarying CDOM and non-algal particles to total light absorption in 
the blue. This large interference from CDOM and non-algal particulate absorption is not 
sufficiently accounted for in the MODIS algorithms. Our results are consistent with findings by 
Harding et al. (2005) who examined the performance of the SeaWiFS empirical chlorophyll 
algorithm OC4v4 in the Chesapeake Bay. The OC4v4 chlorophyll algorithm is the SeaWiFS 
“analog” of the currently default MODIS algorithm OC3M examined here. Wozniak and 
Stramski (2004) used Mie scattering calculations to show that even relatively low concentrations 
of mineral particles of the order of 0.1 g m' 3 can considerably affect chlorophyll estimations from 
standard SeaWiFS and MODIS algorithms that are based on blue-green reflectance ratios. In 
near-shore, estuarine waters like Chesapeake Bay, mineral particles concentrations are typically 
much larger (average concentration of 4.8 g m measured during our cruises; unpublished data, 

G. Gallegos), strongly affecting the Rrs signal in the blue-green. 

Because of the strong interference from CDOM and non-algal particulate absorption in 
the blue-green, regional ly-specific algorithms based on the strong chlorophyll-a fluorescence 
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signal in the near-infrared (NIR) or the chlorophyll-a absorption feature in the red (~ 675 nm) have 
been proposed for improved chlorophyll retrievals in highly turbid, coastal and inland waters 
(e.g. Gower et ah, 1984; Gons, 1999; Ruddick et ah, 2001). When we examined correlation 
between [chl-a] and various MODIS Rts ratios using our in-situ data for the mid-mesohaline 
Chesapeake Bay we found that variability in surface [chl-a] was most strongly correlated with 
changes in the ratio Rrs(677)/Rrs(554) (Table 4). The derived relationship between [chl-a] and 
i?rs(677)/Rrs(554) was consistent with predictions based on RT model simulations and observed 
relationships between in-situ data of [chl-a] and IOPs (b\> and a ) at the two specific wavelengths. 
In the 677 nm wavelength region phytoplankton absorption becomes relatively significant, due to 
the chl-a pigment absorption maximum at ~ 675 nm and the exponential decrease of both CDOM 
and non-algal particulate absorption with increasing wavelength. This, and the almost constant 
red-green ratio of backscattering coefficients measured in the mid-mesohaline Chesapeake Bay, 
suggest that most of the variability in the Rrs ratio 677:554 is due to changes in [chl-a]. 

These results are first step towards the development of regionally-specific chlorophyll 
algorithms for the Chesapeake Bay using remotely sensed ocean reflectance in the red (i.e. ratio 
Rrs(677)/Rrs(554) for MODIS). Moreover, one of the main factors affecting the accuracy of satellite 
chlorophyll retrievals in near shore waters is correction for the atmosphere’s optical characteristics. 
Problems with atmospheric correction of satellite data in the blue due to extrapolation of aerosol 
properties from NIR to shorter wavelengths (Gordon and Voss, 1999) are avoided when using 
the red wavelengths for chlorophyll retrievals. The relationships found here between [chl-a] and 
various Rxs ratios for the Chesapeake Bay strongly suggest that exploiting the ocean color signal 
in the red-NIR, where interference from CDOM and non-algal particulate absorption is minimal, 
is necessary for improving satellite monitoring of biological activity in these Case 2 waters. 
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5.2. Rrs, IOPs, and variability in f/Q 

Since Rrs is related to the backscattering and absorption properties of all the optically 
active water constituents, approximate forms of Eq. 1 have commonly been used in the 1 
development of semi-analytical bio-optical models for remote retrieval of surface [chl-a] and 
IOPs in Case 2 waters. Such approximations are often based on the assumption that the quantity 
f/Q is independent of wavelength, or in-water and air-water boundary conditions (e.g. semi- 
analytic version of chlorophyll algorithm by Carder et al., 2002; Gordon et al., 1988). However, 
recent model results by Morel and Mueller (2002) for vertically homogeneous, Case 1 waters, 
suggest that variations in f/Q (due to its dependence on water IOPs and illumination conditions) 
are within the range 0.08 - 0.15. For nadir-viewing measurements of Rrs, Morel and Mueller 
found that f/Q varies between 0.08 and 0.1 1 depending mainly on [chl-a], wavelength, and solar 
zenith angle (their figure 13.10). 

Our detailed measurements of b b , a, and Rrs spectra allowed us to examine the 
relationship between Rrs and the IOPs ratio b b /(a+b b ), and assess observed variability in the 
quantity fQ for the nadir- viewing geometry of our measurements in the mid-mesohaline 
Chesapeake Bay. Close agreement in the resulting relationships between b b /(a.+b b ) and Rrs, 
whether Rrs was measured in-situ or calculated by the RT model, demonstrated the consistency 
between our in-situ bio-optical dataset and RT model predictions for these Case 2 waters. 

ft t 

Estimated values of the quantity 2 — om) o.»o _ } derived from the linear regression of Rrs meaS ured 

Q n w 2 

versus b b !(a+b b ), ranged from 0.0485 to 0.058. These estimates were in good agreement with the 
ft t 

j — j>. a L (M ‘') values derived based on RT model simulations for the specific bio-optical 
Q nf 
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conditions observed in the Bay waters (i.e. slope of linear regression on i?rs mo d e i versus bd(a+bb). 
Figure 4). 

Since t(a, w /(w,a)/«w is approximately equal to 0.54 (Mobley, 1994), our results regarding 

— • ^"'.a)ha .w) correS p 0n d to flQ values in the range 0.09 to 0.107. Interestingly, the range of jlQ 

Q nf 

values we derived for the Case 2 Chesapeake Bay waters is very similar to the range of f/Q 
values (0.08-0.1 1) reported by Morel and Mueller (2002) in their theoretical computations for 
same geometry of measurements (nadir-viewing), but Case 1 waters. 

Our measurements in the Bay were performed over a wide range of water bio-optical 
characteristics (e.g. [chl-a] in the range 3.5-74 mg m' 3 ), at different solar zenith angles (0 o 
between 16° and 58°), and for both hazy and clear atmospheric conditions. However, coefficients 
of variation in the estimated slopes of the linear regressions between Rrs and b\J{a+bQ (average 
f/Q over a range of conditions) were less than 2.5%, and the standard deviation of the f/Q values 
estimated from individual datasets of Rrs and bdia+bQ was less than 10%. These results, based 
on both in-situ data and RT model simulations, suggest that for the nadir viewing geometry of 
our measurements in these Case 2 waters, observed variability in the quantity f/Q (and therefore 
in the ratio Rrs : bf(a+b\f) was not very large, and not considerably different from Case 1 
waters. 

5.3. Rrs, bb processes, and particulate matter in the mid-mesohaline Chesapeake Bay 

As the spectral reflectance of the ocean is to a first approximation proportional to 
bb/(a+bb), backscattering processes are of primary importance to applications of optical remote 
sensing in oceanography. Backscattering signal depends on the concentration, composition, size 
distribution, shape, and refractive index of suspended, organic and inorganic, marine particles 
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(van de Hulst, 1981). Thus, b\, carries useful information about seawater constituents that affect 
biological activity, biogeochemical cycling, and carbon fluxes in coastal ecosystems. Carbon 
content in individual plankton cells has been found to be related to particle size (e.g. Verity et al., 
1992; Montagnes et al., 1994) and refractive index (Stramski, 1999). Measurements by Stramski 
et al. (1999) in the Southern Ocean revealed high correlation between surface particulate 
backscattering at 510 nm and surface concentration of particulate organic carbon (POC). In 
coastal waters a significant fraction of suspended particulate matter consists of highly refractive, 
inorganic mineral particles, derived from coastal erosion, river discharge of terrigenous inorganic 
particles, bottom resuspension, or aeolian inputs (Stramski et al., 2004). 

Particulate backscattering in the mid-mesohaline Chesapeake Bay showed large variation 
in magnitude, related to variations in particulate loading, mixing processes, and distance off- 
shore. Highly backscattering waters were observed near station JT, which is located closest to the 
land among the four sampling stations and is more strongly influenced by shoreline erosion and 
resuspension of bottom sediments due to tidal and land boundary effects (Figure 2). Variation in 
surface £> b significantly affected Rxs magnitude at all wavelengths (Figure 5), and was found to 
be the main factor driving observed variability in the Rxs signal at 670 nm (Figure 5(c)). The 
strong correlation found between b\> and i?rs(670) indicates that satellite measured Rxs at 670 nm 
can be applied to remotely retrieve particulate backscattering in these Case 2 waters. 

Estimated backscattering fraction in the mid-mesohaline Chesapeake Bay had an average 
value of ~ 0.013 at 530 nm, in agreement with b\,/b values reported by Sydor and Amone, 

(1997) for the near shore waters off Mississippi. The observed spectral shape of bjb is in 
agreement with Mobley et al. (2002) who measured a decrease in b\,/b from 442 to 555 nm by 
less than 24 %, for the Case 2 waters offshore of New Jersey. Backscattering fraction can 
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provide a proxy to the particulate bulk refractive index, which in turn is an indicator of the 
particulate composition in the water (Twardowski et al., 2001; Sullivan et al., 2002, Stramski al., 
2004). Due to their high water content phytoplankton cells have relatively low refractive index 
compared to inorganic particles (e.g. Carder et al., 1974; Aas, 1981 ; 1996; Stramski et al., 1988). 
As particulate backscattering increases with increasing particulate refractive index, b\Jb values 
for phytoplankton-dominated waters are typically lower than those of waters where suspended 
inorganic particles and detrital material dominate (Twardowski et al., 2001; Sullivan et al., 

2002). Measurements by Twardowski et al. (2001) in the Gulf of California and Boss et al. 
(2004) off the New Jersey coast showed that phytoplankton-dominated surface waters with high 
chlorophyll concentrations had bjb values of ~0.005-0.006, while bjb was higher (exceeding 
0.012 according to Boss et al.) in regions where highly refractive re-suspended inorganic 
particles dominated. In agreement with these studies we observed the largest b\Jb values (bjb of 
0.025-0.036) close to the bottom at station JT, consistent with an increase in the relative 
proportion of highly refractive re-suspended inorganic sediments in bottom waters. Surface 
hb/h(530) in the mid-mesohaline Bay had an average value of 0.0125. This is considerably larger 
than the b\Jb values of 0.005-0.006 observed in surface coastal waters dominated by low 
refractive index phytoplankton particles (e.g. Twardowski et al., 2001; Boss et al., 2004). These 
results, and the even higher (average of 0.0146) surface bjb measured at the highly turbid station 
JT, suggest that particulate backscattering in the mid-mesohaline Chesapeake Bay surface waters 
is dominated by suspended non-algal particles, such as highly refractive minerals or organic 
detritus with low water content (Twardowski et al., 2001 ; Green et al., 2003). 

The inference of qualitative geochemical information on non-algal particles from their 
optical characteristics (e.g. absorption magnitude and S nap ) is still a major challenge for effective 
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remote sensing in coastal environments (e.g. Ferrari et al., 2003). The magnitude of non-algal 
particulate absorption in the mid-mesohaline Chesapeake Bay showed strong seasonal and 
temporal variability, and not strong covariation with CDOM or chlorophyll concentration 
(Tzortziou, 2004). However, there was little variation in the spectral shape of the non-algal 
particulate absorption, with S nap values within a narrow range around an average value of 0.01 1 
nm' 1 (SD =0.001 nm' 1 ). Previous studies suggest that S na p values of about 0.01 1 are typical for 
mineral-dominated waters. Bowers et al. (1996) estimated an average S nap value of 0.010 
(± 0.0002) for the absorption spectra of over 100 samples of mineral suspended solids collected 
from the Menai Strait in the Irish Sea. Measurements by Babin et al. 2003) at about 350 stations 
in European coastal waters showed that S nap had an average value of 0.01 17 nm' 1 for the mineral- 
dominated waters in the North Sea and English Channel. Ferrari et al. (2003) measured S nap 
values in the range 0.0095 to 0.0125 for the coastal waters of the North Sea and German Bight 
where the main fraction of total suspended solids was inorganic (~ 76%). Relatively higher S nap 
values, in the range 0.01 15 to 0.0145 were reported for the Baltic Sea, which is known for its 
high organic matter content (Voipo, 1981 ; Ferrari et al., 2003). Similarly, Babin et al. (2003) 
found that highly organic samples collected from the Baltic Sea were characterized by relatively 
higher, and statistically significantly different, S nap values compared to mineral-dominated 
waters, suggesting that observed variation in S nap may be related to the proportion of mineral and 
organic matter. 

The large measured b\>/b values relative to phytoplankton-dominated waters, and the 
close agreement between our estimates of S nap with S nap values reported in other studies for 
suspended inorganic particles, indicate a predominant backscattering role by highly-refractive 
non-algal particles with high inorganic content in the mid-mesohaline Chesapeake Bay. 
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Consistent with these results, we found that b t>, although highly variable, was strongly correlated 
with the magnitude of non-algal particulate absorption at 380 nm (R 2 =0.83, Figure 6). These 
results, along with the small variability observed in S nap , suggest that the high correlation 
between a nap and by, is caused by the dominance of the non-algal inorganic particle 
concentrations in controlling changes in both a nap and by, in these waters. 

Our results of the high correlation between particulate by, and a nap , in conjunction with 
remote retrieval of surface by, from Rrs(670), suggest that Rrs(670) can be applied to remotely 
determine abundance and distribution of non-algal particulate matter in near-shore regions where 
suspended inorganic particles strongly affect ocean color. For the mid-mesohaline Chesapeake 
Bay surface waters we showed that, indeed, J?rs(670) is strongly correlated with the magnitude of 
non-algal particulate absorption at 380 nm (Figure 8). Binding et al. (2003) found good 
correlation between surface irradiance reflectance (R = Eu/Ed) at 665 nm and concentration of 
mineral suspended sediments for the Irish Sea waters. Due to the similarity in the absorption 
spectral shapes of CDOM and non-algal particles, separating contribution by the two 
components to total light absorption in Case 2 waters is a difficult task. As a result, most satellite 
Case 2 algorithms (e.g. Maritorena et al., 2002; Magnuson et al., 2004), including the MODIS 
semi-analytical chlorophyll algorithm (Carder et al., 2002), combine CDOM and non-algal 
particles to one term (e.g. “gelbstoff ’) even though the two components do not covary in coastal 
waters. Regional specific relationships, such as those we have derived here between a nap (380), 
&b(650), and Rrs(670), separately estimate contribution by non-algal particulate matter to total 
light absorption based on the backscattering properties of the non-algal particles and the remote 
retrieval of particulate by, from satellite-measured Rrs. 
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6. Summary and Conclusions 

Effective interpretation of ocean color is a major challenge in near-shore Case 2 waters 
like the Chesapeake Bay. In order to deduce composition and concentration of water 
biogeochemical constituents from ocean color in these waters, it is imperative to apply accurate 
in-situ bio-optical data and RT model results to, first, evaluate the performance of currently 
operational satellite algorithms and, then, develop improved, and in most cases regionally- 
specific, retrieval methods. For the water conditions we observed in the Chesapeake Bay 
currently operational MODIS chlorophyll retrievals are based on empirical relationships between 
[chl-a] and blue-green Rxs ratios. These algorithms cannot be expected to perform accurately in 
Case 2 waters where large concentrations of non-covarying CDOM and non-algal particles 
strongly affect the i?rs signal in the blue. We found that relationships between [chl-a] and two- 
band Rrs ratios in the long visible wavelengths (i.e. i?rs(677)/i?rs(554)), where interference from 
CDOM and non-algal particulate absorption is minimal, provide a better basis for satellite 
monitoring of phytoplankton blooms in the mid-mesohaline Bay. Our measured b\Jb and S nap 
values indicated a dominant role by highly -refractive non-algal particles in regulating 
backscattering variability in these waters, which was consistent with the strong correlation found 
between a nap and particulate b\,. Relating these results to ocean color data, we found that bb is the 
main factor driving observed variability in i?rs(670), a quantity that can be measured remotely. 
Retrieval of particulate bb from satellite measured i?rs(670) can be applied to remotely monitor 
distribution and concentrations of highly refractive suspended particles in the mid-mesohaline 
Chesapeake Bay. Moreover, retrieval of non-algal particulate absorption from the reflectance 
signal at 670 nm could be used in conjunction with inversion of UV-blue wavelengths to derive 
CDOM absorption in these Case 2 waters. Separating contribution by these two similarly 
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absorbing, but non-covarying, components to total light absorption is particularly useful for 
remote sensing of CDOM and studies on dissolved organic carbon cycling in coastal waters. 

Very good optical closure results, achieved using detailed RT modeling, increased confidence on 
the accuracy of, and consistency among the data presented in this study. Quantitative analysis of 
particulate composition, and measurements of the differences in the optical characteristics (i.e. 
Snap and b\>) between organic detrital and inorganic mineral suspended particles are needed to 
improve interpretation of remote sensing in coastal regions like Chesapeake Bay where non-algal 
particles often dominate the backscattered signal. 
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Figure Legends: 

Figure 1: Location of in-situ measurements (stations HB, PI, TI and JT) and cruise track. The 
starting point during our cruises was the SERC dock located in the Rhode River sub-estuary, 
along the western shore of the Chesapeake Bay. 

Figure 2: Frequency histogram of measured surface [chi -a]. 

Figure 3: Comparisons between measured [chl-a] and chlorophyll concentrations derived from 
the MODIS empirical algorithms using our field data and estimates of Rrs and nZw spectra 
as input to the algorithms (chlorophyll values in logarithmic scales; N=40). The MODIS chl- 
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algorithms used were: (a) “Chlor_a” (Eq. 4); ‘chlor_a_3’ (Eq. 6) (b) for the ‘unpackaged’, 
and (c) for the ‘fully packaged’ cases; (d) ‘chlor MODIS’ (Eq. 5) for the high chl-a pigment 
(always the case in mid- and upper-Chesapeake Bay); (e) ‘chlor MODIS’, 5th order 
polynomial (D. Clark, updated 19 March 2003)). 

Figure 4: Relation between measured Rrs(k) values and measured — Mjj — (solid circles), and 

a (X)+b b (X) 

between model-estimated Rrs(k) values and measured — — (open circles) at (a) 440 

a W+b„(X) 

nm, (b) 530 nm and (c) 670 nm. Linear least-squares regression fits are shown as solid lines 
for i?rs measured and dashed lines for i?rs model (R 2 between 0.90-0.98; Pd.0001 in all cases). 
Regression results are summarized in Table 4. 

Figure 5: Relation between measured b b and measured Rvs (solid circles), and between measured 
b b and model-estimated i?rs (open circles) at wavelengths (a) 443 nm (b b measured at 450 
nm), (b) 532 nm (b b measured at 530 nm) and (c) 670 nm (b b measured at 650 nm) (R 2 are 
given in Table 4). 

Figure 6: Relation between surface measurements of 6 b (650) and a nap at (a) 380 nm (R =0.83, 
N=43, P<0.0001) and (b) 412 nm (R 2 =0.76, N=43, P<0.0001). 

Figure 7: Relation between surface [chl-a] and surface 6 b (650) measurements in the Chesapeake 
Bay waters (R are given in Table 4). 

Figure 8: Linear regression between measured or model-estimated (based on measured IOPs) 

Rrs values at 670 nm and surface measurements of absorption by non-algal particles at (a) 

380 nm and (b) 412 nm (N=33, P<0.0001) (R 2 are given in Table 4). 
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Table 1: Dates of cruises in the Bay, and instrumentation for measurements of radiation fields. 
The AC9 and ECOVSF3 instruments were used in all cruises for measurements of water IOPs 


Date of cruise Instrument used 

for radiation fields 


2001, June 4 
2001, June 11 
2001, June 25 
2001, July 9 
2001, September 21 
2001, September 26 
2001, September 28 
2001, October 4 
2001, October 30 

2001, November 13 

2002, May 6 
2002, May 15 
2002, May 22 
2002, June 6 
2002, June 18 
2002, June 28 
2002, November 8 


Satlantic OCI-200 
Satlantic OCI-200 
Satlantic OCI-200 
Satlantic OCI-200 
Satlantic SMSR 
Satlantic MicroPro 
Satlantic MicroPro 
Satlantic OCI-200 
Satlantic MicroPro 
Satlantic OCI-200 
Satlantic MicroPro 
Satlantic MicroPro 
Satlantic MicroPro 
Satlantic OCI-200 
Satlantic OCI-200 

Satlantic MicroPro 
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Table 2: In situ measurements of IOPs and radiation fields in the Chesapeake Bay. 
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"Remote Sensing Reflectance and Inherent Optical Properties in the Chesapeake Bay " 


Table 3: Summary of M0D1S empirical chl-a algorithms 

MODIS empirical chl-a algorithms 

Parameters 

MODIS OC3M, “Chlor a” product: 

logio[chl-a] = A+ BR 3M + C-R 3M 2 + DR 3M 3 +E-R 3M 4 

(Eq 4) 

fl3M=logio(max[ * rs(443) , ^ r5(488) ]) 

Rrs( 551) flrs(551) 

A=0.283 
B=-2.753 
C=1.457 
D=0.659 
E=- 1.403 

MODIS empirical-HPLC, “chlor_MODIS” product: 

log 10 [chl-a]=[A(log 1 o(X)) 3 +B(log 1 o(X)) 2 +C(log 1 o(X)) 1 +D]/E 

(Eq 5) 

Different sets of A-E parameters are used for waters with high and 
low [chl-a]. nLw is the normalized water leaving radiance (Gordon 
and Clark; 1981), which can be estimated by multiplying 7?rs(A,) with 
the extraterrestrial solar irradiance, Fo(7,)). 

v _nLw(443)+nLw(488) 

nLw(551) 

High Chl-a pigm Low Chl-a pigm. 
A=-2.8237 A=-8.1067 

B=4.7122 B=12.0707 

C=-3.9110 C=-6.0171 

D=0.8904 D=0.8791 

E=1.0 E=1.0 

Carder et al (2002), “chlor_a_3” product: 

logio[chl-a]= c 0 + C]logio(X)+c 2 (logio(X)) 2 +c 3 (logio(X)) 3 

(Eq 6) 

Parameters c 0 - c 3 are adjusted dynamically (based on information on 
the sea surface temperature) in order to account for pigment packaging 
effects (‘packaged’, ‘un-packaged’, and transitional cases) in nutrient- 
replete and nutrient-deplete conditions. 

Rrs{ 488) 
i?w(551) 

Unpackaged Fully Packaged 
c 0 = 0.2818 c 0 = 0.51 

ci= -2.783 cj= -2.34 

c 2 = 1.863 c 2 = 0.4 

c 3 = -2.387 c 3 = 0.0 
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"Remote Sensing Reflectance and Inherent Optical Properties in the Chesapealce Bay " 


Table 4: Regression Relationships between IOPs and Rrs in the mid-mesohaline Chesapeake 
Bay 


Parameters (Y vs X) 

Regression Relationship 

R 2 

jRrs(412)//frs(554) vs [chl-a]j„. S i, u 

logio(Y) = -0.354 • logio(X) - 0.4508 

0.24 

7?rs(443)/7?rs(554) vs [chl-a] in . situ 

logio(Y) = -0.1537 ■ logio(X) - 0.2933 

0.39 

Rrs(488)//?rs(554) vs [chl-a], n . situ 

logio(Y) = -0.1299 • logio(X)- 0.1074 

0.32 

i?rs(510)/Rrs(554) vs [chl-a] in . S itu 

logio(Y) = -0.1 168 • logio(X) - 0.0339 

0.36 

7?rs(532)/JRrs(554) vs [chl-a] i„. S it u 

logio(Y) = -0.0728 • log 10 (X)+ 0.0068 

0.23 

7?rs(670)//?rs(554) vs [chl-a] , n . S i tu 

logio(Y) = 0.166 - logio(X) - 0.5467 

0.48 

i?rs(677)//?rs(554) vs [chl-a] m . sim 

log 10 (Y) = 0.1725 ■ log 10 (X)- 0.5117 

0.54 

^t-w(^77) VS [chl-a] i n -situ 

Y = 0.0166 -X + 0.0603 

0.92 

a t . w (554) vs a t . w ( 677) 

Y = 0.68 -X + 0.08 

0.92 

b b (6 50) vs 2>„(530) 

Y = 0.74 ■ X 

0.99 

i?rs(443) vs ^(450) 

Y = 0.0496 ■ X, for i?rs(443) mode i 

0.94 

a w (443) + a,_ w (443) + b b (450) 

Y = 0.0485 ■ X, for Rrs(443) data 

0.90 

/?rs(532) vs ^ (530) 

Y = 0.0533 ■ X, for Rrs(530) mode i 

0.96 

a w (530) + a,_ w (530) + b h (530) 

Y = 0.0498 • X, for flrs(530) data 

0.98 

i?rs(670) vs ^(650) 

Y = 0.057 ■ X, for 7?rs(670) model 

0.96 

a w (670) + a,_ w (670) + b b (650) 

Y = 0.058 • X, for Rrs(670) data 

0.94 

6 b (450) vs Rrs(443) 

Y = 34.06 • X - 0.02 

0.67 

b h ( 530) vs Rrs(532) 

Y = 12.63 -X -0.017 

0.77 

h b (650) vs Rrs(670) 

Y = 15.82 -X -0.008 

0.88 

Z> b (650) vs a nap (412) 

Y = 0.0463 • X - 0.0092 

0.76 

&b(650) VS «nap(380) 

Y =0.0411 ■ X- 0.0107 

0.83 

[chl-a] vs fe b (650) 

Y = 433.5 -X -0.92 

0.57 

a nap (412) vs 7?rs(670) 

Y = 302.2 -X + 0.159 

0.70 

a nap (380) vs 7frs(670) 

Y = 356.71 -X + 0.168 

0.74 
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Figure 4, Tzortziou et al.. 
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Figure 5, Tzortziou et al.. 
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“ Remote Sensing Reflectance and Inherent Optical Properties in the Chesapeake Bay 


Popular Summary: 

The particulate and chemical composition of the ocean and their concentrations can. be 
measured from satellites by monitoring the color of the water. Interpretation of the measured 
color is still a major problem,, especially in coastal waters such as the Chesapeake Bay. Recently, 
we have shown how to combine theoretical models for the propagation of light through water 
with an extensive bio-optical database to examine the relative roles of light absorption and 
scattering from different substances that might be in the Chesapeake Bay. We applied the results 
to measurements just above the water’s surface and to satellite using our derived relationships 
between remotely sensed quantities and optical properties of individual water components. We 
also evaluated methods for remote retrieval of the amount of chlorophyll and suspended 
inorganic particles. 

Currently operational chlorophyll algorithms used by MODIS (MODerate resolution 
Imaging Spectroradiometer) are based on empirical relationships between chlorophyll 
concentration ([chl-a]) and blue-green remote sensing reflectance (Rrs) ratios. These algorithms 
break down in coastal waters because they do not sufficiently account for the effects that 
substances other than chlorophyll (i.e. non-living particles and dissolved organics) have on the 
color of the water. We show that relationships between [chl-a] and Rrs in the long visible 
wavelengths, where interference from non-living particles and dissolved organics is minimal, 
provide a better basis for satellite monitoring of phytoplankton blooms in these waters. Both our 
measurements and model simulations indicate that highly-refractive non-living, mainly 
inorganic, particles play a dominant role in regulating backscattering variability in the 
Chesapeake Bay waters. These results, in conjunction with our findings that backscattering in 
these waters can be estimated from the satellite measured Rrs at 670 nm, suggest that satellite 
ocean color can be applied to remotely monitor distribution and concentrations of highly 
refractive inorganic suspended particles in the Chesapeake Bay. Also, our results axe the first 
showing that satellite measurements can be used to separate the contribution by non-living 
. particles and dissolved organics to the total light absorption in the blue wavelengths. Separating 
the contribution by these two separate, but similarly absorbing, components to total light 
absorption is particularly useful for studies- on dissolved organic carbon cycling in the coastal 
zone. 
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